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Abstract We prove the existence of static solutions to the cylindrically symmetric 

Einstein- Vlasov system, and we show that the matter cylinder has finite extension in 

L^ two of the three spatial dimensions. The same results are also proved for a quite general 

class of equations of state for perfect fluids coupled to the Einstein equations, extending 

the class of equations of state considered in [3]. We also obtain this result for the 

H , Vlasov-Poisson system. 



1 Introduction 

The most frequently studied static spacetimes containing matter are the spherically sym- 
metric ones. This class of spacetimes is compatible with asymptotic flatness. Another 
interesting class of asymptotically flat spacetimes are the axially symmetric ones. How- 
ever, a study of this class is at present out of reach. As a first possible step towards axial 
symmetry it could be useful to study the cylindrically symmetric spacetimes, although 
cylindrical symmetry is not compatible with asymptotic flatness, nor are there any cylin- 
drical configurations in the universe. Nevertheless, cylindrically symmetric spacetimes 
can be considered as a subclass of the axially symmetric ones with an extra symmetry 



requirement. It should also be pointed out that time dependent cylindrically symmetric 
spacetimes admit gravitational waves, in contrast to the spherically symmetric case, and 
still reduce the complexity of the Einstein equations. 
^ In 1917 Levi-Civita obtained the most general static cylindrically symmetric vacuum 

solutions. Since then many explicit cylindrical solutions have been obtained in the fluid 
case, cf. [8]. Unfortunately these solutions are often local, and no global analysis is 
usually available. However, in [3] the global properties of static cylindrically symmetric 
spacetimes with perfect fluid matter is studied, and global existence of solutions and 
finiteness of the radius of the fluid cylinder is shown. It is shown that when the fluid 
cylinder has finite extension, it is possible to glue it smoothly with a Levi-Civita solution, 
and in this way obtain a global solution. To prove the finite extension of the fluid cylinder 
the authors make use of the fact that the boundary density, i.e. the energy density of 
the fluid on the boundary where the pressure vanishes, is positive in an essential way. 
The main part of this paper is concerned with the Einstein- Vlasov system. However, we 



will also improve previous results for perfect fluids by including all equations of state 
that satisfy 

< °o, (1) 



/o Pip) + P 

where p is the energy density, p the pressure and P c the central pressure. The 
Einstein- Vlasov system is the Einstein equations coupled to a collisionless kinetic matter 
model, cf. [TJ E] for an introduction to the Einstein- Vlasov system. Prior to this study, 
only the spherically symmetric Einstein- Vlasov system has been studied in the static 
case, cf. [6j [13] HH] I15j . However, in the time dependent situation more symmetry 
classes have been studied, cf. [T] for a review of the results on the Einstein- Vlasov 
system. Cylindrically symmetric time dependent solutions have been considered in [5]. 

For the static spherically symmetric Einstein- Vlasov system the proof of local exis- 
tence is straightforward, but the proof of global existence is more involved. The proof 
of finiteness of extension of matter is essentially done in one of the following ways. Ei- 
ther it is proved by using the non-relativistic limit, knowing that this limit has finite 
extension of matter, or it is proved by introducing new variables in an appropriate way, 
and by showing that these functions blow up at a finite radius one concludes that the 
matter has finite extension, cf. [13] and [15] for the different approaches. We will rely 
on the last approach. We point out that in [B] a completely different approach is used to 
investigate the static spherically symmetric Einstein- Vlasov system. In [6] the system 
of equations is transformed to an autonomous dynamical system on a state space with 
compact closure, and the solution orbits are studied to extract information whether the 
corresponding solutions have finite radius or not. 

For the static cylindrically symmetric Einstein- Vlasov system the proof of local ex- 
istence is straightforward, but uses the fact that the matter terms are continuously 
differentiable. To ensure this differentiability property and to be able to compare our re- 
sults with the spherically symmetric case, and also to be able to handle the matter terms 
more easily, we study the following class of distribution functions. Let f(r, v) := 4>{E)L ■ , 
I > — 1 , where E is the particle energy and L is the modulus of the angular momentum, 
cf. section 2. We also require that 4> £ ^((0, Eq)), where Eq > is a constant, a cut-off 
energy. To obtain global existence we need to do a careful analysis of the equations 
involved. The proof of finite extension of matter do not need the asymptotic structure 
of the distribution function, in contrast to the spherically symmetric case, where it is 
needed. Heuristically, one of the reasons that the proof of finiteness of extension is eas- 
ier in the cylindrically symmetric case than in the spherically symmetric, is that the 
singularity of the equations at r = is stronger in the latter. 

As mentioned above we also prove the same result for the static cylindrically sym- 
metric Einstein equations coupled to a perfect fluid which satisfies equation ([1]). In 
the analogous results for the spherically symmetric case the asymptotically polytropic 
behavior as the pressure tends to zero is taken into account, which is not needed for 
the cylindrically symmetric case. This is in turn analogous to the comparison between 
the cylindrically- and spherically symmetric Einstein- Vlasov system as mentioned above. 



We also show that the matter solutions, both for Vlasov matter and perfect fluids, can 
be glued with a unique vacuum solution. 

In the final part of the paper we also consider the non-relativistic Vlasov- Poisson sys- 
tem. Global existence have been proved previously in [2j under more general conditions 
than we have, but we show in addition that the matter cylinder has finite extension. 

2 Static cylindrical spacetimes 

Consider a spacetime with two hypersurface orthogonal Killing vectors, where one is a 
translation £ z , and the other a rotation £#. Let the spacetime also be equipped with a 
timelike Killing vector ^. Then the most general metric can be written in the form, cf. 



ds 2 = _ e 2u dt 2 + e 2( 7 -^ r 2 + e 2^2 + ^-2^ 

where t £ (—00,00), r £ [0, 00), z £ (—00,00) and 9 £ [0, 2tt), and where u, 7 and 
ip only depend on r. We normalize the gravitational constant G = 1, and the speed of 
light c = 1. Then the Einstein matter equations read, cf. [16J, 

t = ^* _ 8vre 2 ^-^T/, (2) 

^±^ =V ,'2 +87re 2( 7 -^) r;) (3) 

r 

ip" + — + ifj'87rre 2 ^-^ (T r r + T l t ) = Aite 2 ^-^ {T e e - T z z + T r r + if) , (4) 

v" + - + ^Svrre 2 ^-^) (T r r + T/) = 4vre 2(7 - ,/ ' ) (T r r + T e 9 + T* - T\) , (5) 

where T?, i = t,r,8 and z are the components of the energy-momentum tensor. For 
a perfect fluid the nonzero components of the energy-momentum tensor read 

Tt = - P , r; = Tl = T e e = P, 

where p is the energy density and P is the pressure. From V a T fe a = we obtain 

^ = -v>( P + P). (6) 



In kinetic theory the energy-momentum tensor is defined by 

/ T) T)h 1 

T ab (x) = - f(x,p)-?—y/\g\dp, 

Jr3 po 



where a, b = 0, 1,2,3 and \g\ is the determinant of the spacetime metric, and x are 
the spatial coordinates with corresponding momentum p, i.e. p is the spatial part of the 



four-momentum. Furthermore p° is expressed through p a and x via g a bP a P b = — 1> where 
we have assumed that all particles have the same rest mass m = 1. We introduce new 
momentum variables v a , a = 0, 1, 2, 3 by 



v ° = e v p°, v 1 = e 7 "V, v 2 = eV, v 3 = re"V, 



and note that 



v 







y/i + (v 1 ) 2 + (v 2 ) 2 + (v 3 y 



The nonzero components of the energy-momentum tensor for static cylindrically 
symmetric kinetic matter then read 

T\ = - f f(r,v)v°dv:=-p, 

Jm.3 

Tl= \ f(r,v) { -^dv := P u i = 1,2,3. 

We immediately obtain the following important inequality 

p>P 1 + P 2 + P 3 , (7) 

which holds on the r support of /. The distribution function / satisfies the Vlasov 
equation 

v'df , Av 2 )*,., lJv 3 ) 2 df 
,v 1 v 2 df 1 , v 1 v 3 df 

-V<^r^-(--VO^r^3=o. (8) 

Now multiply the Vlasov equation (jHJ) by v 1 and integrate in velocity space, assuming 
compact support in the velocity variables, or at least sufficiently fast decay to zero as 
\v\ — ► oo, then we obtain 

HP 1 

( L± = - u \ p + Pl )-l l l,'P l + -P l . ( 9 ) 

dr r 

Note that this reduces to equation (JBJ) in the particular case of a locally isotropic 
distribution function, i.e. the distribution function only depends on the energy, see 
below, since then P\ = P2 = P3 ='■ P- 

We have three constants of motion. The particle energy E = e u v°, the angular 
momentum squared L 2 = r 2 e~ 2 ^\v 3 \ 2 and Z = e^v 2 . If f(r,v) = $>(E,L,Z), with 
& *E C 1 , then / satisfies the Vlasov equation. We will however restrict the function class 
and only consider functions of the type 

f(r,v) = cf ) (E)L l , 1>-1. 

Furthermore, <fi G C 1 ((0, Eq)), Eq > and k > —1 such that on every compact subset 
K C (0,co) there exists a constant C > such that < </>(E) < C(E - E)\. The 



matter terms will be continuously differentiable under the assumption k + -^- > 0, cf. 
Lemma 1. Then the matter terms read 

P(r) = ^ m £ HE)E*(E*-e^dE, (10) 

rEo 
Pi(r) = ATTA ir l e^ l+A)v ^ / (f>(E)(E 2 - e 2v )*dE, i = 1,2,3, (11) 

Je v 

where 

1 

A x = A 2 



(Z + l)(Z + 3)' 
and 



A, 



1 + 3 
Now define 

"-Bo 



H m (u):= / cj>{E){E 2 -u 2 ) m dE, 

J u 



then (1101) and dill) become 



4 7rr « P -('+4)^-^ 

I + 1 2 2 

Pi(r) = 4vr^r'e-^ +4 )^ /,/, J H' 1 +3(e 1/ ), i = 1,2,3. (13) 

2 

The boundary conditions which supplement equations ©-© are given by 

dl dv , dib „ . ,. 

7 = -i = !/ = — =j, = -Z = 0, at r = 0, 14 

ar ar ar 

and follow from the regularity on the axis r = 0, cf. [3]. 

We will also consider the non-relativistic counterpart to the Einstein- Vlasov system, 
namely the Vlasov-Poisson system which in the static cylindrically symmetric case reads 

-(rU'{r))' = 4vrp, (15) 

r 

Here p is given by p(r) := J" K3 f(r,v)dv. The conserved quantities in the non- 
relativistic setting are the energy E := 2M 2 + U{r), the angular momentum squared 
1? := r 2 |u| 2 sin 2 #, where u := (f 1 ,^ 2 ) and 6 the angle between x and u in the plane 
of fixed z, and the quantity Z := v 3 . As in the relativistic case we have that any rea- 
sonable function which only depends on these variables satisfies the Vlasov equation. 



By restricting to the same class of functions as before, i.e. f(r,v) := <j)(E)L l , I > — 1, 
4> G C l {{— oo, Eq)), Eq > — oo a constant, and 4> satisfying the requirements in Lemma 
1 below, we obtain 

1+5 

p( r ) = ^- rl 9l±l(U(r)), (17) 

l + 1 2 

where 

/•-Bo 
<7m(«):= / 4>{E){E-u) m dE. (18) 

Define £ m j n := -co in the non-relativistic case and E m i n := in the relativistic case. 
The following lemma will be needed later, for a proof cf. [15J. 

Lemma 1. /73]/ Let <j> : (E m i n , oo) — > R 6e measurable, Eq > E m i n , and k > —1 sitc/i 
t/mi on every compact subset K C (E m i n , oo) i/iere exists C > sitc/i i/iai < ^>(-E/) < 
C(£'o — E)\,E G ET. Lei H m (u) and g m {u) be defined as above i.e. 

pEo 

H m (u) = I 4>{E)(E 2 - u 2 ) m dE, u G (0, oo), 

Ju 

g m (u) := / cp(E)(E - u) m dE, u G (-oo, oo). 

Ju 

If m > —1 and k + m + 1 > then g m ,H m G C((E m i n , oo)). Ifm>0 and k + m > 
then g m ,H m G C l {(E min) oo)) w#j 

5m = -mg m -i, 

H' m = -2muH m -i. 

With a distribution function <fi as in Lemma 1, we have that when u > Eq, then 
H m {u) = g m ( u ) = 0, which imply that the matter terms vanish. 

3 Existence and finite extension of matter 

3.1 The Einstein- Vlasov system 

For later purposes it is convenient to rewrite the Einstein equations with Vlasov matter 
as follows. We follow [16] and introduce the new variable 



M :-- 



I(l_ e 2("+tf-y)). (19) 



The reason for this transformation is that the system then has a structure similar to 
the static spherically symmetric Einstein- Vlasov system, with M(r) used in a similar 



way as the local ADM mass m(r) in the spherically symmetric case, cf. |15j . Note that 
cylindrical spacetimes cannot be asymptotically flat and hence do not have any ADM 
mass. The Einstein matter equations are transformed to 

— = 2nre 2v {p-P x ), (20) 

dib 47T [ r 9l . 1 . „ . , , . 

' re 2v ■= = (P 3 -p)dr, (21) 



dr rVl - 8M Jo VI - 8M 

? = /" r re 2 - - * (P 3 + 2P X + p)dr, (22) 

dr rVl - 8M 7 Vl-8M V r/ V ; 

where we have used the regularity at the axis. 

The boundary conditions follow from equation fl!4f) and read 

,, dM dv , dib 

M = —— = u = — =ip = -j- = 0, at r = 0. 23 

or ar ar 

We have the following local (in r) existence theorem. 

Theorem 1. Let f(r,v) = <j){E)L ,1 > —1 and 4> as in Lemma 1 with k + -^- > so 
that H m G C l ((0,oo)), m = -^-, -^. Then for 5 > sufficiently small, there exists a 
unique C 2 -solution, (pf,ib,v), to the system (p|)-([5j), (Q with boundary conditions (\14h 
on the interval [0, 5] . 

Theorem 1 follows immediately from the next lemma, since the system ©-©, © 
with boundary conditions (J14H is equivalent to system (I20p -(l22"j). ([8]) with boundary 
conditions (1231). 



Lemma 2. Lei / and </> 6e as above, then for 5 > sufficiently small, there exists a 
unique C 1 -solution (M,ip,u) to the system (fl?0j) — (f^gjl on £/ie interval [0,5]. 

Proof. Take 5 > 0. Define the set 

B = {w = (M,ib,v) G C°([0,o],]R 3 ) :M(0) =0,-0(0) =0,i/(0) = 0, 

|| M ||oo< Ki, || V Hoo< K 2 , || 1/ ||oo< ^ 3 }, 

where K{,i = 1,2,3 are positive constants, K\ < |, and B is equipped with the 
norm || w ||b:=|| M ||oo + || t/> lloo + || v \ 00 . B is a closed set in the Banach space 
C ([0, o],IR 3 ), and thus a complete metric space. 

Now define the operator T, acting on B by 

Tw:= (fTw)i,(Tw)2,(Tw) 3 ), 

where 



rr 

(Tw)i := 2vr / fe 2v {p - Pi)df, 
Jo 



(24) 



(iy 



2 ■' 



4tt 



(Tw] 



3 • = 



o rVl - 8M Jo 

4vr 

re 



re 



2v 



2,' 



1 



V1-8M 
1 



(P3 — p)dfdf, 



P 3 + 2Pi + p)^^. 



(25) 



(26) 



/o fVl - 8M 7o V1-8M 

Recall that the matter terms p, Pi,P2,P 3 are functionals of the metric components 

and hence has a w dependence. First we show that T(B) C B. Take w £ B. Since 

w G P, we have that e^ S [e - ^ 3 ^^ 3 ]. Recall by Lemma 1, Hi+i,Hi+3 G C 1 ((0,cxd)), 

2 2 

and since the arguments take values in a closed bounded interval which do not contain 

zero, H1+1 ,Hi+3, H' l+1 and H' l+3 are all bounded. The argument is similar for all the 

2 2 — — 

components of the mapping, so we only show that the second component of the mapping 

satisfies the claim above. Since I > — 1, equations (fT2j) and (fl"3j) yields the estimate 



re 



2// 



1 



VI - 8M 



(P 3 -p)|<C7r 



HI 



This implies that 



4tt 



~~VT 



f e 2u 

8M Jo VT 



8M 



(P 3 - /9)dr 



< C7f' +1 , 



J+2 



which yields 

II (Tu») 2 ||oo< Cr l 

Similar estimates yield that || (Tw)i ||b< Cr' +2 , i = 1,2,3. Hence, for 5 small enough 
Tw £ B, and since w was arbitrary T(B) C B. We point out that the constant C above 
may change from line to line and does not depend on 5. 

We now claim that T acts as a contraction on the space B for 5 sufficiently small. 
Then, by the Banach contraction theorem we have a unique continuous solution to the 
system (I20p - (]23p . but it is then straightforward to verify that this solution in fact is a 
C 1 -solution. We want to estimate || Tw — Tw \\b where w,w £ B. Since the estimates 
of the different terms are similar we just sketch the proof of the second term. The mean 
value theorem will be used to estimate the differences, and Lemma 1 will be used to 
ensure the continuity of the matter terms and its derivatives. Now the difference in the 
second component of the mapping T can be written 



\(Tw) 2 - (Tw) 2 \ 



-in 



fy/1 - SM 



re 



2;/ 



1 



VI - 8M 



(P3 — p)dfdf 



4tt 



re 



20 



v\J\ - 8M Jo Vl-8M 



:(p3 — p)dfdf 



4ir 



f V1-8M v/i_8M 



re 



:2v 



1 



+ 



4tt 



f\J\- 8M Jo 



f[(e 



2v , e 2P, 



1 



V1-8M 



+ e"( 



V1-8M 
1 1 



V1-8M Vl - 8M 



(P 3 - p)df 

)](P 3 -p)dr 



h VI - 8Af 



(P 3 - P 3 - p + p)df] 



df 



where P3 and p indicates the dependence on w. Define 

F{v,Tjj):=e-tt l+2 >+W, 
then, since w,w 6 B, by the mean value theorem 



FM)-F{v,$)\ < sup (e-W+y°+Wy/(l + 2) 2 + l 2 }\(v,il,)-(i>,4,)\ 



u,ip€B 



Now define 



< C || w — w \\b ■ 
G{M) :-- ' 



VI -8M' 
and again by the mean value theorem 

|G(M) - G(M)| < C || M - M ||oo 

< C || u> — u; ||b • 
By using equations (fT2J) and (fT3l) we obtain as before that 

1 



Ire 



V1-8M 



(P3-P) <Cf 



d+i 



Hence, 



r 4vr 1 



1 = )/V 1 



/o f V1-8M Vl-8M Jo V1-8M 

< Cr' +2 || 10 -u> ||s . 
Similarly by using the mean value theorem we have 

I 1v 2u I ■ /~i 11 ~ 11 

e — e I < C -|| to — w ||s, 
which yields 

4vr f f . 1 



(P3 — p)dfdf 



fVl - 8M Jo \/l-8M 

< (V +2 || w — to 



(P 3 - p)((? v - e 2v )drdf 



and 



, 47r I fe 2i> (P 3 - p)( . 1 - , 1 )dfdf 

f VI - 8M Jo VI - 8M Vi _ 8 m 



(27) 



(28) 



< Cr l+2 || w - w \\ B ■ (29) 
By the mean value theorem again 

\H m (en-H m (e )\<sup\H' m (en\\e»-e°\ 

< C || w — w ||s, 
with m = -4p or m = -jp, which implies that 

/ , . / re 2i> - „ (P 3 ~ ^3 ~ P + P)drdf 

./o fVl-8M-/o V1-8M 

< Cr'+ 2 || u; - w \\ B ■ (30) 
Equations (|27j) - ([30]) yield 

|| (Tw) 2 - (Tw) 2 ||oo< Cr^ 2 || «; - u) || fl . 

Similar estimates hold for the first and third component of the mapping T. Hence 
for 5 > small enough, T is a contraction mapping on B since 1 + 2 > 1. Thus, there 
exists a unique continuous solution to the system (|20p — (|22p on the interval [0,5], with 
^(0) = i/(0) = M(0) = 0. Since I > -1 we obtain by using the equations pD>p2j) that 
the solution in fact is a C 1 -solution that verifies the boundary conditions (|23p . 

□ 

Remark. By using equations ©-© we immediately obtain that (y, tp, 7) is a C 2 -solution. 

Define 

p ■= ^e-^+^-^Hi+s^), 

2 

then Pi = AiP, i = 1,2, 3. The following lemma gives us some useful information about 
the behavior of the matter and the metric. 

Lemma 3. Suppose we have a C 1 -solution (M, ijj, u) to the system ( \2(M - ft23§ with f(r, v) = 
4>{E)L l , I > — 1, k + £t- > and 4> o,s in Lemma 1, then M , v and 7 are strictly in- 
creasing and ip is strictly decreasing as long as p > and M < | . 

Proof. From equations (|20l) . ([2]) and (|22l) we immediately obtain that M, 7 and v are 
strictly increasing since p > 0. Furthermore from equation (|2ip we obtain that ip is 
strictly decreasing. D 



Define 
and 



R er := maximal interval of existence, 



R := sup{r < R ex : p(f) > 0, for f £ [0, r)}. 



10 



Lemma 4. Assume that there exists ro G (0,R ex ) such that p(ro) = 0, then p(r) = 
forr G [r ,R ex ). 

Proof. If /o(ro) = then e u ^ r °' > Eq and by monotonicity, e u ^ r > > Eq and hence /)(r) = 
for r > r$. □ 

It is important to control M as the next two lemmas show. 

Lemma 5. If r £ [0,R ex ) then M{r) < |. 

Proof. Recall the definition of M which reads 

On the domain of existence v, ip and 7 have finite values, and hence M < g. □ 

We conclude that if M(f) = h then the solution breaks down at r = f. The next 
corollary follows immediately. 

Corollary 1. Let f(r,v) := cf)(E)L l , I > -1 and k + ^±± > 0. If p{R) = for some R, 
then a unique vacuum solution, can be joined to the matter cylinder. 

Proof. We have a unique solution on [0,12], since R < R ex . Since equations (|20]) - (|22]) 
are valid in the vacuum case too, and since all functions and their derivatives have finite 
limits as r — > R~ , we can take these values as boundary values for a unique Levi-Civita 
solution. □ 

If R ex = +00 the solution is global (in r), which in fact is true as the next theorem 
shows. 

Theorem 2. The system (TJ2j)-(rj3j), with I > — 1, k + ^ > and <j) as in Lemma 1, 
has a unique global (in r) solution, i.e. R ex = +00. 

Proof. Assume that R ex < +00. First we show that R ex = R. Assume that 12 < R ex , 
then we can glue a unique vacuum solution to the matter solution, but then the solution 
can not blow up at finite R ex as assumed. Hence, 12 = R ex . 

Now we want to show that M E C 1 ([0,12]). First, e v < Eq since v is monotonically 
increasing. By Lemma 3, ip and v are monotonically decreasing respectively increasing, 
so by equation ([3]) and the boundary conditions (|14l) we see that |^| < v. Hence, tp and v 
have finite limits as r -> R~ , and ^f = Ivre 2 " ^*'^^ [H1+3 (e v ) + e 2u H l±1 (e u )] - 
AnA^e-^^- 1 ^ Hi+3(e u ) has a finite limit as r -> R~ , since H m E C 1 , m = ^-, ^. 

Hence, M G C X ($,R\), and we can expand l-8M(r) = (R-r)S(r), with S G C([0,12]). 
Next we prove that M(r) < C < |, r G [f ,12). By equations ([H]) and ([22]) we 
observe that \ip'\ and 1/ are bounded if M < C < 1/8. Hence by equation ([2]) 7' is 
also bounded and the solution can be extended contradicting the assumption that 12 ex 
is the maximal interval of existence. Recall by Lemma 3, Lemma 4, and M(0) = 0, that 



11 



M(r) /O, r£ [ijf,R), since otherwise there is no matter in spacetime and it becomes a 
Levi-Civita spacetime with regular axis, i.e. Minkowski spacetime. Now define 

4tt f^^-PiK- 4tt r A*lr^e-^>-^Hi^{e») 

CJ "= / — -(IT = / - df 

ry/l - 8M Jo V1-8M " ry/l - 8M Jo (/ + 1)(Z + 3)V1 -8M 

(31) 
then equations (J2TJ) and (J22j) can be written as 

2 (M + 4vrr 2 e 2 "Pi) + Q[rQ - -=j ]. (33) 



dr r(l-8M) v " ^ L ^ v/1 - 8M 

We have to consider two cases. 

D QTirl rMmno ( ) I vl I — _ . , , , , , , , , — 



First let -1< I < 0. Then Q < and hence Q[rQ - , x ] > 0. With r£[§,R) 



-(M + 47rrV' y P 1 ) + Q[rQ- ]> 2 - 



dr r(l-8M) v 7 ^ L ^ VI - 8M ~ 1 - SM ~ (R - r)S(r) ' 

and integration yields that lim u(r) = +oo which contradicts the fact that v(r) is 

uniformly bounded when r G [0,-R). Hence, lim M(r) < g, and the solution can be 

T-+R- 
extended, contradicting the fact that R ex is the maximal interval of existence. Thus the 

solution is global (in r) when — 1 < I < 0. 

Now consider the case / > and define 

7o Vi - 8M 

By equation ([3]) and the boundary conditions (|14l) 

fe 2l/ Pi 

rrjrfr < i/(r) + ^(r) < oo, r G [0, i?). 

q i — oivi 

Since P3 — P2 = IP\ we have 

Q := f ^W-^j _ , f _j£S_* < 00, r £ [0, R). 

Jo V1-8M 7o V1-8M L 7 

We observe that Q(0) = and Q'(r) > 0, r G (0, i?). 
Now by equations ([3]) and (l32j) 

^> r ^-^ = l[l- J_^ ] 2 + rQ 2 + Q[l-^=L=]- : l[l-^=L=]-Q 
ar 4r Vl - 8M Vl - 8M 2r Vl - 8M 

' 'i - 4vrQ + 16tt 2 Q 2 ) - -L. 



r(l-8M)M ' 4r 

12 



Thus if Q(r) < C < ±, r G [0, R) we have that \ + 16vr 2 Q 2 - 4vrQ > C > 0. Hence 

— > - - — (34) 

dr ~ r{R — r)S(r) 4r' 

and integration of this inequality from % to r and letting r tend to R will lead to a 
contradiction. 

To obtain this, observe that the integrand in equation (I2ip . has the same sign so 
ij)' < 0, r € [0, R) , and equation (j32j) can be written as 

^ = — + (-1 + 8ttQ), r£ (Q,R). (35) 

2r 2rVl - 8M 

Since Q is increasing, lim Q(r) = C. Now there are three possibilities. First of all 

r^R~ 

assume C < g^, then we are done since this is actually what we want to prove. 

Secondly assume that C > g^. Since Q is continuous, this means that Q(r) = ^, 
for some r < R. However, then by equation (I35D ip'(r) > 0, which is a contradiction. 

Hence, we only have to exclude the case C = lim Q{r) = i. Assume by contradic- 



»il- 



tion that lim Q(r) = g^, and consider 



-1 + 8ttQ . . 

(36) 



V1-8M 



If lim , _ =f > — 1, then ij)'(f) > for some f G [0, R) by continuity, which is a 



contradiction. Now by PHospital's rule 



.. -l + 8vrQ ,. 8vrQ' 
lim — . = lim 



,. v/l-8A/ ,. "l - "3 
= Imi — - — 5—7 — =rr = nm -— 

r ^R- -tore*" (p-F\) r ^ R - p - Pl 

V1-8A/ 

-IH^ie") _j 

= rSS- (Z + 2)ifw(e v ) + (/ + 3)e^H 1±1 ( e ") " (T+2) > ~ L 

2 2 

Hence we obtained a contradiction, which implies that 

Q(r)<C<^, r €[*,*). 

Now integrate the inequality (I34D from | to r, and let r tend to i? we see that v 
diverges to +00, as r — > 1?~, which again is a contradiction. Hence M(r) < C < |, r G 
[0, i?) and again the solution is global (in r). 

□ 
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The following corollary is an immediate consequence of Theorem 2. 

Corollary 2. Let f(r,v) := 4>{E)L l , I > —1 and k + l -^- > 0. Then there exists a 
unique global (in r) C 2 -solution to the system (0) — ((5j) and Q subject to the boundary 
conditions < \14\) - 

With the class of distribution functions considered above we could have two possibil- 
ities. Namely, either the matter fills all the space or there is a finite radius of extension. 
In fact as the next theorem shows, the first possibility above is excluded for the class of 
distribution functions considered in this paper. 

Theorem 3. Let f(r,v) := (j)(E)L l , I > —1 and k + ^±1 > 0. Then there is a finite 
radius R such that p(r) = 0, r > R. 

Proof. Assume by contradiction that R = +oo. Since we have matter (otherwise R = 0), 
and the matter terms are continuous, we can assume that there exists t\ £ (0, +oo) such 
that p{r\) 7^ 0. By continuity of the matter terms p(r) ^ 0, r £ {r\ — e, n + e) for some 
e > 0. By equation (I20p and Lemma 4 



M(r) > e, (37) 

for some e > when r >r\. We also have 

1 - Vl - 8M > C > 0, r G [n, +oo). (38) 

Define 

7 ? :=(logS -^)+, (39) 

and 

«==^. (40) 

V 

Thus ry is a well defined and differentiable function since logE^ := vr, see equation (JX2J) . 
Furthermore, 

rf = ~y, (41) 

where ' denotes differentiation with respect to r. Since n(0) = vr > 0, £ is also a well 
defined and differentiable function when r £ [0,-R). Now 

, M' Mrf M' Mv' Mv' v'~ 

« v ( ry )2 v ^ (^2 - ( ?? )2 M « • ^ ^ 

By equations ([^U|) and (f2"2"j) we have 

dl/ 47T /' r _ 9 „ 1 



dr ry/1 - SM J y/l- SM 



[ f ^7fhM iP ^ 2Pi+p)df (43) 



vra i (1 " i/r3W) + 4I [ ^^ha^ + 3n)df] (44) 
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> 1-yrag (45) 

~ 2rVl - 8M 
By equation (|3gj) . 1- ^/l - 8M(r) > C > 0, r G [n,R), and thus 

c > -e. (46) 

r 

Integrate this inequality from r\ to r. Then we immediately obtain that there exists 
R < +oo, such that £ — » oo as r — > P~. But this implies that 77 tends to zero since M{r) 
is bounded and increasing. Hence, there exists an R < 00 such that p(r) = 0, r > R 
which contradicts the assumption above. □ 

3.2 Perfect fluid 

In this section we will consider the Einstein equations coupled to a perfect fluid. By 
definition, a barotropic fluid is one whose energy density and pressure are related by an 
equation of state that does not contain the temperature. We will assume the fluid to 
have a barotropic equation of state 

P = P(P), (47) 

where p and P are the energy density respectively the isotropic pressure of the fluid. 
We will also assume that the equation of state p = p(P), is C 1 when P > 0, and 
p(P) — > pr < +00, P — > 0, where pr > is a constant. 

Many of the results in this section is already known, cf. |3j. However, for com- 
pleteness we sketch the proof of the results. The proofs essentially follows the proofs 
in section 3.1 for the corresponding results. The important new result of this section is 
that we are able to treat a more general class of equations of state than what is done 
in [3]. Hence Theorem 5 is the main result of this section. The following local existence 
result follows from Lemma 7 below. 

Proposition 1. Let p(P), be a barotropic equation of state, p £ C 1 ((0,P c ]) ; where P c is 
the central pressure and furthermore limp(P) = pR, where pr, < +00 is a nonnegative 

constant. For each value of the central density p c = p(P c ) > 0, there exists r$ > 0, such 
that a unique C 2 -solution of the field equations d2|)-([6]) with boundary conditions (\14h 
exists for r £ [0, r<j] • 

In the perfect fluid case with M := g(l — e 2 '^" 1 "^ -7 )) as in section 3.1, the equations 
corresponding to equations (f20|) - ([22|) can be obtained by observing that the perfect fluid 
case corresponds to the I = case for Vlasov matter. Hence, P\ = P2 = P3 := P in this 
case, and by using equation (|20p it is possible to integrate equations (|2ip and (|22p to 
obtain equations (fiU]) and (foTH) . The analogue of equations (f2"U|) - (f2"3"|) with equation © 
added are 

— = 2Trre 2 »(p-P), (48) 
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d ± = Lu- _ 1 1 (49) 

dr 2r [i VT^8M 1, [ ' 

dv 2 ,.. . 2 2u 



dr r(l - 8M) 



(M + 47rrV"P), (50) 



W = ~^ P + p) ' (51) 

M = ^ = , = ^=V = ^ = 0, atr = 0,P(0)=P c >0. (52) 

dr dr dr 

Then we have the following result. 

Lemma 6. Let p(P), be a barotropic equation of state, p G C 1 ((0,P c ]), where P c is 
the central pressure and furthermore lim p{P) = Pr, where pr < +oo is a nonnegative 

constant. For each value of the central density p c = p(P c ) > 0, there exists r^ > ; 
such that a unique ^-solution of the equations ([7^|)- (f57]) with boundary conditions (f5l|) . 
exists for r £ [0, ro] . 

Proof. The proof follows the proof of Lemma 2 with / = 0. □ 

From equations ([2|), and (|48p — (|5ip we can easily obtain the following lemma. 



Lemma 7. Suppose we have a solution with a regular axis. Then M , v and 7 are strictly 
increasing and ip and P are strictly decreasing as long as P > 0. 

Proof. Using equations ([2]), (|48p -(|5ip we can immediately read off the claimed properties 
of the different functions. □ 

The following is an easy consequence of equation (|51|) . We recall the definitions 

R ex := the maximal interval of existence, (53) 

and 

R := sup{r < R ex : p{f) > 0, for f G [0,r)}. (54) 

The following lemma holds also in the perfect fluid case. 

Lemma 8. Assume that there exists r^ £ (0,R ex ) such that p(ro) = 0, then p(r) = 
forr £ [r ,R ex ). 

Proof. Let ro be as above. Then p(ro) = and since P < p, P(ro) = 0. By equation 
(f5i~T) . P is a decreasing function, thus P(r) = 0, r G [ro,R ex ). But p considered as a 
function of P implies that p(r) = 0, r G [ro, R ex )- D 

Lemma 5 and Lemma 6 holds without any changes in the proofs. 

Lemma 9. If r G [0,R ex ), then M(r) < |. 
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The analogue of Theorem 2 holds for perfect fluid matter since we have equation 
(f50j) which reads 

At/ 9 _ _ 

-(M + 4vrr 2 e 2iy P). 



dr r(l - 8M) 

Then it is easy to apply the same strategy as in the proof of Theorem 2. We have the 
following corollary which says that we can extend the spacetime by adding a unique 
vacuum solution to it. 

Corollary 3. If P(R) = 0, for some R < oo ; then a unique vacuum solution can be 
joined to the fluid cylinder. 

Proof. See the sketch of proof of Corollary 2. □ 

Theorem 4. Let p{P), be a barotropic equation of state, p £ C 1 ((0,P c )), where P c is 
the central pressure and furthermore 1hap(P) = pn, where pn < +oo is a nonnegative 

constant. Then for each value of the central density p c = p(P c ) > 0, there exists a unique 
global (in r), C 1 -solution of the equations (\48\)- l\51h with boundary conditions (j5^j) . 

The following corollary is an immediate consequence of Theorem 4. 

Corollary 4. Let p(P) be a barotropic equation of state, then there exists a unique global 
(in r) C 2 -solution to the system (fjj) — ([TJ) subject to the boundary conditions fli^p and 
< P c < +oo. 

Now we have two possibilities. Either the fluid fills all the space or there is a finite 
radius where the fluid vanishes. 

In [3] they proved that for perfect fluids with positive boundary densities p(0) > 0, 
and regular axis, there always exists R < oo, such that P(R) = 0. This fact can be 
generalized to all equations of state which satisfy 

dp 

< oo. 



/o Pip) + P 
Hence, we have the following theorem 

Theorem 5. For a solution with regular axis and an equation of state that satisfies 
po = p(P c ) > and 

' '' < oo, (55) 



/o pip) + V 
there is a finite radius R such that P(R) = 0. 

Proof. The proof follows the proof of Theorem 2, but we outline it for completeness. 
Define 

Jo pip) + P 
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and 

n 

Observe that £ is well defined, since it is matter at the axis of symmetry and the equation 
of state is continuous. Furthermore rj is well-defined and r/(P c ) is finite and nonzero. We 
also have that 77 is continuous, and M is finite. By (15 ip 

7/ = -V, (58) 

where ' denotes differentiation with respect to r. Now M(r) > e, r G [n,+oo), e > 
(see the proof of Theorem 3, section 3.1) so by equation ((501) 

77 (77) 2 77 (r?) 2 ~~ [j]) 2 M ~ r 

From this inequality we immediately obtain that there exists R < 00, such that £ — ► 00, 
as r — > i?~. But this implies that r\ tends to zero since M < g. Hence, there exists an 
R<oo, such that P(R) = 0. D 

Remark. Observe that 77 is not bounded for linear equations of state of the type P = 

(p—l)p, 1 < \i < 2 which is in complete agreement with the explicit solutions obtained in 

p i_ 1 

[7]. In [16] however they consider polytropic equations of state p = (-^) M +nP, n = ^~— 

where K and p are constants, and these models have finite radius as can be seen by 

evaluation of the integral (|56p in the definition of 77. Also observe that for equations of 

state with positive boundary density, i.e. p(0) > 0, the integral ((56|) is finite, compare 

with the results in [3]. 

Remark. The proof of finiteness of the radius is much easier than in [TO] because we only 
need to consider one equation and all terms turn out to have the right sign. However, 
if one tries to apply exactly the same strategy in the spherically symmetric case it fails 
due to an r _2 -term that appears in the differentiation of rj. 

3.3 The Vlasov-Poisson system 

Consider now the system (j!5j) — (I16p . which with the ansatz f(r,v) = 4>{E)L < , I > —1, 
and 4> satisfies the assumptions in Lemma 1, is turned into a nonlinear equation for U 

1+9 o 

(rU')' = 2 -^-r l+l gi±1 {U), (60) 

t + 1 2 

where 

g m {U):= / <f>(E)(E-U) m . 

Ju 

Let U satisfy the boundary conditions U(0) = a, U'(0) = 0, where a is a constant. 
Then we have the following local existence result, cf. [2] where they prove a more general 
existence theorem. 
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Proposition 2. Let f(r,v) = (fi(E)L 1 , I > —1, and <f> as in Lemma 1 with k + ^fci > 
such that gi+i 6 C 1 (( — oo, +oo)). Then there exists 5 > 0, such that equation (EEI), 

2 

subject to the boundary conditions £7(0) = a, U'(0) = 0, has a unique C 1 -solution when 
re [0,5]. 

Proof. The proof is omitted since it is similar to the proof of Lemma 2. □ 

The solution is global as the next theorem shows, cf. |2j for a more general theorem. 

Proposition 3. Let f(r,v) = <j)(E)L l , I > —1, and 4> o,s in Lemma 1 with k + ^jp > 
such that gi+i G C 1 ((— oo, oo)). Then there exists a unique C 1 -solution of equation (fgQj) . 

2 

subject to the boundary conditions U(0) = a, U'(0) = 0, on [0, oo). 

Proof. Let R ex and R be defined as in the relativistic case. Assume that R ex < +oo. 
Clearly U is increasing. Either r — ► R ex U{r) < E$ on its maximal interval of existence, 

lim 

which implies that 

gi+i(U)<C, (61) 

2 

by Lemma 1, and by equation (|60f) 

U'(r) < Cr l+1 , (62) 

which is uniformly bounded on [0, R ex ), since I > —1. Thus the solution can be extended 

which contradicts the maximality of the solution, and hence the solution is global (in r). 

If on the other hand r — ► R ex U(r) > Eq, r £ [0, R ex ), then for all r > R, p(r) = 0, 

lim 

and again U exists globally since we can add a vacuum solution to the matter solution 
at r = R. □ 

In fact with / as above the matter has finite extension which we show below. 

Theorem 6. Let f(r,v) := <j){E)L , I > —1, k + -^- > and cj) be as in Lemma 1, then 
there exists an R, such that p(r) = 0, r > R. 

Proof. Define 

pr 

m := 2tt / fp(r)df, 
Jo 

and 

n := (E - U)+, 

then m and r\ are well-defined and satisfy the following equations 

m = 2irrp, 
and 

n 



: 2m 



r 
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where ' denotes differentiation with respect to r. Now define 

m 

4 : = — , 
V 

then £ satisfies 

m' mrf 2 2 

^ = 7^1 - ~£ • 

7/ (r/) z r 

As in the relativistic case, there exist an R such that £(r) — ► oo, as r — ► I? - . Hence 
there exists i? such that p(r) = 0, r > R. D 

Acknowledgment: I want to thank Hakan Andreasson for useful discussions and 
comments. 
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